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1. Introduction 


Thomas-Fermi theory or semiclassics is а convenient method to formulate average 
ground state properties of nuclei, like matter and level densities, ground state energy 
or kinetic energy density; for an overview, see Ring and Schuck[1]. The same methods 
were also applied to properties of excited nuclei, like multiphonon level densities(2,3], re- 
sponse functions[4], or the optical potential and quantities derived therefrom[5]. Only a 
few applications, however, treat good angular momentum like Jensen and Luttinger[6| 
who calculate angular momentum distributions in the Thomas-Fermi model. Other 
quantities like multiphonon excitation level densities and response functions are also 
experimentally accessible for given £ and correlated occupation numbers and the ef- 
fective mass of f-states could be compared with those of a Brueckner-Hartree-Fock 
calculation. 

In this paper, we therefore formulate semiclassics for given angular momentum and 
given z-component or summed over the z-components. In section 2, we derive the 
propagator for given £ т and other quantities. In section 3 the transformation to 
Wigner space is discussed. Section 4 contains applications to one-particle and one- 
particle one-hole level densities and to the imaginary part of the optical potential and 
the last section summarizes the paper. More details of the theory can be found in ref.[7]. 
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2. Semiclassics in Coordinate Space 


In semiclassical approximation the eigenfunctions of the Hamiltonian locally are plane 


waves 
<= >= е" | (1) 
and 
Ak? у 2 
< НН >= —— + V(r) (2) 


are the eigenvalues of the Hamiltonian. For spherically symmetric potentials which 
commute with @ and @, the projection onto good angular momentum / and its z- 
component m is achieved by employing the eigenstates of Р and i, 


< r|k£m >= 4ri“ j(kr) Ур (ДУ, (Р), (3) 
or, summed over the magnetic quantum numbers, 
< rikl >= (2£ + 1)i! (Кт) Pr(cos), (4) 


where созӣ = К. r/kr. Here and in the following indices Ё, m or £ denote fixed mod- 


ulus and z-component of the angular momentum or summed over the z-component, 
respectively. In semiclassics, the single particle propagator in coordinate representation 


C*(r, г!) =< rleF# |r’ > (5) 


is the central quantity. Here, Û is a Lagrange parameter to constrain the energy, which 
mostly serves as Laplace-transformed energy, ёз. For В = 0 one obtains the unit 
operator 

Cp, r') = 6(r — r’). (6) 


For fixed £, m or £ this becomes 


C3>°(r, x) = Pim (r, r)ó(r ~r’), 
Се, т) = P(r, r')é(r — т). (7) 


Here we have defined the projectors 
Pr, r) = Te mec) 


_ 26+ 1 P{cos x) 


Pis. v 
(rr) = 217089 (8) 
where cos x = r-r'/rr’. The same result is obtained from the constraint 

< r|(2£ + 1)6(2 — (Е + 1)) 660, — m)|r' >= Pem(r, x’) 6(r — "7: (9) 


As a result, the semiclassical traces become 
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T =4 [ае fer б(т—т), 
Trim = 4 Је" far Pan (r, F') 6(г-г'). (10) 


The factor of four comes from summation over spin and isospin whereas in tbe local 
quantities (not under traces) this is not taken into account. In treating the potential as 
locally constant the propagators can be evaluated by neglecting h-corrections, i.e. using 
the asymptotic form of the Bessel function, see appendix, with the argument mr?/8^?, 


Ск) = (25) е (ame 2 via) , 


1/2 CH 
Zmrr 


а 0 
?кћ?р / 


{om f. \ 
С? (кг -( Pin (к, к’) exp [- —35(r-r) -8giv(R)-- n). 11 
ter.) ter ¥) exp [gt did. سف‎ )) (11) 
Here and in the following 44+1) is replaced by / in the spirit of lowest order semiclassics 
but 204-1 is left because it facilitates summation over £. From the propagator one derives 


the nonlocal spectral density 


g'(r, 1) =< r|ó(e — Н)|т' >= 251, C" (r), (12) 
the partition function 
29 = Tr C?(r,r), (13) 
the density of states e 
gle} = Tr 6(e — Н) = £51, 2°, (14) 
and the nonlocal density 
рвут!) =< re = B)! >= £54, C" (1/8, (15) 


where Л is the Fermi energy measured from the bottom of the potential. Here, relative 
and center-of mass coordinates are defined as 
rtr 

2 LJ 
and the potential У (Е) is normalized to V(R =0) = 0. With the propagators (11) the 
nonlocal densities becomes 


В = 


s-r-r, (16) 


яе} = 8) dh Cru , 


i — r|pc (R)/^ 
Pte (x, T’) = 1р0) e con) д (17) 


where the local Fermi momenta are defined as 


pr(R) = y2m(A – V(R)) Ө(А – V(R)), 


РЕК) = 2m (a-ve) e (à-vm- 55). 


(18) 
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and the same local assumption has been made as in the potential, ie. rr' has been 
replaced by R?. From this one derives the local densities 


them А (19) 


which vanish beyond the classical turning point R, where pr(R,) = 0 and the kinetic 
energy densities 


_ Folr,r') РК) 
В) = эгэр lemen = рзд 
(4 3 
Tem( R) = = 0 (20) 


Since the quantities given above do not depend on the magnetic quantum number т, 
the corresponding m-summed quantities simply obtain by multiplying by a factor of 
(24 + 1). The single-particle level densities reads 


E a [ar e V(r) e - vin, 


9'() = St “Ут f dr (-v- eye (-v- 2) у (21) 


The limits of integration are determined by the classically allowed region defined by 
the step function. The level density for given angular momentum can also be found in 


ref.[6] in integrated form. 


3. Semiclassics in Wigner space 


Semiclassics is also conveniently formulated in Wigner space (R, P), where the relative 


coordinate s is Fourier transformed into P, 


JR, P) = [з “Fh f(x, a). (22) 
The trace becomes is 
Tr - 4 ? E KF (23) 


and the local limit is reached by integrating over all momenta. Since the Wigner trans- 
form of the Hamiltonian is equal to the classical Hamiltonian, the relevant quantities 


obtain simple forms. The propagator, for instance becomes 


С(Е, Р) = eae, (24) 
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where H(P, R) = P*/2m + V(R) is the classical Hamiltonian, and the nonlocal density 

reads 
/(R, P) = © (à - H(P, R)). (25) 


Quantum mechanically and also semiclassically only £ and т are good quantum num- 
bers whereas classically all three components of the angular momentum are fixed. This 
means that the quantum average over the z, y-components of the angular momentum 
is represented by {һе azimuthal average 


$0 -1)-4 / Ф, L) 5(1— T)  &(£— £)5(m — т). (26) 


The classical limit then obtains from[8] ¢,,4, < €, = €. The Wigner transform of eq.(9), 
hence is given by the classical counterpart 


&(R x P/h — 1) = ô(RP sin0/ħ — £) 6((R x Р), /h — т), (27) 


where cos? = R- P/RP. Then the trace becomes 


Tro. = 4 fer PR F(R x P/h — 1). (28) 


T ar 
Since in Wigner space and neglecting h-corrections the expectation value of the product 
of two operators is equal to the product of the expectation values, all single particle 
quantities become proportional to 6(R x P/ħ— 1). The propagator for instance be- 


comes 
С? (R, P) = $(R x P/h — 1) e-?F(^41, (29) 
which, summed over m, can be rewritten as 
pi BA? \ 
СИВ,Р) = 5(РЕзш@/\ — £) exp -8 (3 SP uA RA vu) , (30) 


This now contains the familiar radial kinetic energy by virtue of Рр = Р cos6 and the 
centrifugal potential. The nonlocal spectral density 


Jim(R, P) = $(R x P/A — 1) 6(e - H(P, R)) (31) 
simply reflects energy and momentum conservation and the nonlocal density reads 
ры(В,Р) = 8(R x P/A – 1) O(A — H(P, R)). (32) 


Experimentally, however, only angle averaged quantities are accessible. Hence, averag- 
ing over the polar angles 0, ¢ the angular momentum constraint becomes 


EEIE j= М O(RP/ħ— 4) 


4 2RP\/R?P2/h? — 2 


The nonlocal density, for instance, ede over all angles, becomes 


6(R, Р) = (33) 
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< eR, P) >= 6'(R, P) Ө(А - P*/2m — V(R)), (34) 
from which one derives eq.(19) by momentum integration. 


The quantities derived above by Wigner transforming eq.(9) can also be obtained 
by direct Wigner transform of the semiclassically asymptotic quantities (see appendix). 


4. Densities and level densities 


In this section, we calculate densities and level densities and the imaginary part of 
the optical potential for given angular momentum by employing an harmonic oscillator 
potential. The oscillator constant is denoted by Aw = 41MeV/N'/5, and R, = 1.2£m N'/5 
is used for the radius. The Fermi energy follows from number conservation, 


А = (3N/2) ^ hu. (35) 


The partition function can be evaluated to give 


Z = 151 «е (36) 
and the level density becomes 
2641 
= ө(5-0). (87) 


The level density of the harmonic oscillator simply consists of step functions for given 
energy. The local densities of an harmonic oscillator for fixed angular momentum is 
shown in Fig. 1. 
One notices that only the £ = 0 states contribute at the center. The maximum angular 
momentum for given energy follow from the step function in eq.(18), 
Le) = e/hw. (38) 
Next, we calculate the one-particle one-hole level density. Integrated over all angular 


momenta it is defined as 
Фр (€) = Tr, Tr; O(A — H;) Ө(Н, — A) 6(e — Н, + H2). (39) 
For fixed £, m this goes over into 
di (€). = Tr Tr; f dé, fdl, f dm, [dm,8(A — H+) €(H,—2) &(e— H+ He) 
x6(R, x P,/À — 1 5(R; x P2/h — 12)6(1— (lı + h)). (40) 


By symmetry, the result does not depend on m and the dependence of the azimuthal 
angles of the traces is solely contained in the last delta function. Eq.(40) is therefore 


rewritten as 
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Figure 1: 
Spin-isospin summed density of 200 particles for given angular momentum (4 = 0...6). 


dia (c) = Tr, Tr ав, f dhf dm, Јат.Ө(А- H;) Ө(Н,-А) 6(e— H, + H2) 
х6{ В, Ри sin 0, /№ — £,)6(R; P, sin 0,/ћ — t) А(416,, £), (41) 
where 
a44, 4= | sqq, cya 
_ 10244-12 - 6 - 6) T 
* (466 - (P - 6 - Gp 
is the angle averaged triangular relation. 
By using the methods of refs.[2,3] it is transformed into a folding expression for fixed 
angular momentum 
t bin ГА m 
Hale) = fde, [de fdt, [at в) єє – е +) (6,0. (43) 
Here 9i» and ф{ are the one-particle and one-hole level densities, i.e. the single particle 


level densities restricted to energies above and below the Fermi energy, respectively, 
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Figure 2: 
One-particle one-hole excitation level density of 200 particles for given angular momen- 


tum (£ — 5...20). 


vine) = 9'(«) Ole — А), 
gin) = g'() Ө(А - е). (44) 
Note that our result eq.(43) differs from the expression recently proposed by Кгарре[9] 
who did not .account for the angle averaging in the angular momentum conserving 
delta function but it is conform with the semiclassical limit of the response function of 
Delafiore[10]. The result is shown in Fig. 2. The levelling off which starts at the Fermi 
energy and the later onset of higher angular momenta is due to the maximal angular 
momenta at the Fermi energy. 
Finally, the on-shell value of the imaginary part of the optical potential including 
polarization and correlation graphs is calculated with a zero range two-body interaction 


u(r —r') = —vjó(r — r). (45) 
Without angular momentum constraints in coordinate space it is obtained from the 
imaginary part of the mass орегаќог[5], 
W(e rr) = $12 £5. (« ripe |! >< ripe ^ir >? 
+ < вре #Н|ү' >< r|pe?* |r’ >?) , (46) 
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where 2, are the particle and hole density operators. In semiclassical approximation 
for fixed £, m in Wigner space it becomes 
w^" (e, R, P) = v 29. fap, fdp [Фр,Ө(р}(В)—Р})Ө(р1—Р(В))Ө(р}—Р(Е)) 
(2) 
x&(2m(e—A) - pe (R)--pi - Fi -P3)6 (px -- P2 P3 — P)A(R x (pı pa p3)/A —1) (47) 


Since the force (45) conserves the momentum and, furthermore, since it puts all intrinsic 


coordinates at the same position R it conserves also the angular momentum. Аз a result, 
Ит (e, В, P) remains a single particle quantity and, hence, 


(є, R, P) = W (+, В, P)5(R x P/A — 1). (48) 
The on shell value of the latter has been calculated in ref.[5] as 
0 1-09 < 2me < р - Ph 
» — Е E ‚РЕ = ра < 2me < ph 
v = Ce Pals T п аа at © 
ch. TR ,P PR < 2me < оо 


where C = v;m/(2xh^)* and pg = pp(R) is the local Fermi energy of eq. (18) and the 
on-shell condition is Р2/2т = є - V(R). With the angle averaged angular momentum 
constraint the imaginary part of the optical potential for fixed £ becomes 


W'(e, R) = &“( Po, R} W (e, R). (50) 


As strength for the two-body interaction, v = 745 MeV fm? is used which reproduces 
the experimentally known quadratic behaviour around the Fermi energy and a harmonic 
oscillator potential is employed. The result is shown in Fig.3. 


One observes that due to the classical condition 
ve 
2mR? 
for energies below the Fermi energy only a few angular momenta contribute. The higher 
partial bound states have a larger Wt and since the single-particle spreading width is 
proportional to the imaginary part of the optical potential also a larger spreading width. 


e> V(R) (51) 


This seems to be in agreement with recent observations of Langevin-Joliot et. al.[11] 
At high energy the imaginary part of the optical potential behaves like the centrifugal 
potential. 
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Figure 3: 


Imaginary part of the optical potential of 200 particles at R=5 fm in a harmonic oscil- 
lator for given angular momentum (£ = 2...16). The Fermi energy is at 47 MeV. 


5. Summary and outlook — 


Semiclassics has been generalized to cases where the the modulus and the z-component 
of the angular momentum is fixed. We derived the propagator and other quantities 
in (r,r’) - space and in Wigner space (R,P). In coordinate space quantities under 
traces are easily evaluated because angular momentum projection can be done solely in 
the traces. In Wigner space all single-particle quantities are proportional to an angle 
averaged angular momentum conserving delta function. Ав applications, densities and 
one-particle one-hole excitation level densities have been calculated for an harmonic 
oscillator potential and the imaginary part of the optical potential has been obtained 
for a zero range two-body interaction and with an harmonic oscillator potential. The 
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schemes outlined are very flexible and can be applied to arbitrary potentials and work 
is under way to employ realistic Woods-Saxon potentials and also deformed potentials 
which do not commute with LS 

Furthermore, higher multiphonon excitation level densities will be calculated by 
using the analytic methods of refs.[2,3] as well as response functions. Also the real 
part of the optical potential for fixed £ and the therefrom derived correlated occupation 
numbers and the effective mass are of great interest in order to calculate realistic nuclear 
charge densities({12]. 
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Appendix: Mathematical details 


In the appendix we give some relevant formulas for the derivation of the asymptotic 
propagators. 
Eq. (11) is derived by using the asymptotic expansion of the modified Веззе! function 
of large argument 
Ленц) — (212) ^ е4, (52) 


which semiclassically is the first term of an A- expansion in the large quantity 1/№. 


Direct transforms from coordinate to Wigner space of the £- projected quantities 
proceed as follows. First, the asymptotic expansion of the Legendre polynomials for 
large £ in the projector of eq.(8) 


PF (cosa) — C7 С" „(+ 1/2)x) (53) 


is used to transform the projector from a space to (R,s) space, 


РЕ, в) = ir xa Jo 5s sino), (54) 


where созо = В.-в/Ёз. In taking the R-axis as integration axis it is noted that in 
cylindrical coordinates в = (s, = s sing, з, = s cosa, ¢), the projector depends only on 
55, whereas the remaining integrands depend only on s,, for instance 


5(r — r’) = &(s,). (55) 


Then the Fourier transform (s — P) separates and the integral representation of the 
Bessel function 
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' n 
Jats) = f see сенін (о-и /2) “ 
2 


е 
and the orthogonality relation of the Bessel function, 


Tas s Jy(us)J,(vs) = = Й (57) 


are used to reduce the Wigner transform to a one-dimensional Fourier transform, 


[seh (zs. fs.) = 6(PRsino/n=2) fds, ёте" f(s,). (58) 
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STRESZCZENIE 


Sformuzowano wyrazenia semiklasyczne dla ustalonego modutu 
1 z-owej składowej momentu pędu. Jako funkcji falowych użyto, 
wyrzutowanych na określony moment pędu, fal płaskich, uzyskując 
semiklasyczne wyrażenia dla jednocząstkowego propagatora i wy- 
nikajacych z niego wielkości, jak: partycja funkcji, nielokalna 
macierz gęstości, gęstość poziomów jednocząstkowych i gęstość 
pozioméy czastka-dziura. Y przestrzeni Wignera ва one proporcjo- 
nalne do kombinacji momentu рейџ: 68 x Р\/в 209) A 

ER x $),/n = ш), Formelizm ten zastosowano do wyliczenia 

cZastka-dziura oraz urojonej części potencjału optycznego na 
bazie oscylatora hamonicmego z 2-ciarowym oddziaływaniem o za- 
siggu zerowym. 


РЕЗЮМЕ 


Сформулированы полуклассические выражения для определенно- 
го значения и 2-составляющей момента количества движения. В Ka- 
честве волновых функций применены плоские волны проектированные 
на заданный момент количества движения. Получены полуклассичес- 
кие выражения для одночастичного пропагатора и вытекающих из не- 
го величин: распределение функции, нелокальная матрица плотно- 
сти, плотность одночастичных и частично-дырочных уровней. В про- 
стоанстве Вигнера они пропорциональны комбинации момента коли- 
чества движения: 5(1& xP|/n-1) и Rx P), Ё -ш). Этот форма- 
лизм применялся для расчета плотности одночастичных уровней и 
состояний частица-дырка, а также ынимой части оптического потен- 
циала на основе гармонического осциллятора с двухтельным взаимо- 
действием с нулевым радиусом действия. 
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